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ABSTRACT 


Program  SCATSPHERE2,  a  modification  of  program 
SCATSPHERE  by  F.M.  Henderson,  computes  the  series  solution  for 
the  time-harmonic  plane-wave  scattering  by  a  submerged,  elastic, 
spherical,  evacuated  or  fluid-filled  thin  shell.  For  program 
SCATSPHERE2,  this  report  presents  the  basic  scattering  theory,  the 
user’s  instructions,  and  three  examples  involving  underwater  plane- 
wave  scattering  by  (1)  an  evacuated  steel  spherical  thin  shell,  (2)  a 
fluid-filled  steel  spherical  thin  shell,  and  (3)  a  fluid-filled  aluminum 
spherical  thin  snell. 


ADMINISTRATIVE  INFORMATION 

The  work  was  sponsored  by  the  DARPA  structural  acoustics  project  at  the  David  Taylor  Research 
Center  funded  under  DARPA  Funding  Document  GD32i0890T06660,  Amendment  11,  Program  Element 
63569E.  The  DTRC  program  manager  was  J.R.  Peoples,  Jr  (Code  1930.6).  The  DARPA  program  manager 
was  Dr.  A.  Tucker. 


INTRODUCTION 

A  few  years  ago,  the  NASHUA  structural  acoustics  procedure  1-5  was  created  to  model  the 
radiation/scattering  problem  for  submerged  elastic  obstacles  by  using  the  boundary  element  method  for  the 
exterior  fluid  and  the  finite  element  method  for  the  structure.  The  NASHUA  procedure  could  also  model 
interior  fluids  using  structural  finite  elements  by  using  analogies  between  the  equations  for  dynamic  elasticity 
and  acoustics.  6-7  To  verify  the  NASHUA  procedure,  three  benchmark  problems,  where  the  analytical  series 
solutions  8-11  are  known,  were  chosen:  (1)  time-harmonic  uniformly-driven  radiation  by  a  submerged  elastic 
spherical  thin  shell;  (2)  time-harmonic  sector-driven  radiation  by  a  submerged  elastic  spherical  thin  shell;  and 
(3)  time-harmonic  plane-wave  scattering  by  a  submerged  elastic  spherical  thin  shell.  Programs 
RADSPHERE  12  and  SCATSPHERE  were  developed  by  F.M.  Henderson  to  compute  the  series  solutions 
for  the  second  and  third  benchmark  problems. 
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Recently,  the  NASHUA  procedure  was  enhanced  13  to  model  structural  acoustic  problems  for 
submerged  fluid-filled  structures  by  using  the  boundary  element  method  for  both  the  exterior  and  interior 
fluids,  and  the  finite  element  method  for  the  structure.  For  interior  fluids,  using  the  boundary  element 
method  instead  of  the  finite  element  method  reduces  the  matrix  computation  significantly  for  fine  meshes.  To 
verify  the  new  procedure,  the  time-harmonic  plane-wave  scattering  by  a  submerged  elastic  fluid-filled  spherical 
thin  shell  was  chosen  as  a  benchmark  problem  since  the  analytical  series  solution  can  be  derived  and 
computed.  Here,  a  thin  shell  assumption  is  used  to  simplify  the  analytical  series  solution.  Program 
SCATSPHERE2,  a  modification  of  program  SCATSPIIERE,  was  developed  to  compute  the  series  solution 
for  scattering  by  a  submerged  elastic  spherical  thin  shell  that  is  either  evacuated  or  fluid-filled. 

In  the  next  section,  the  time-independent  scattering  theory  is  introduced.  A  brief  review  of  the  plane- 
wave  scattering  by  a  submerged  evacuated  spherical  thin  shell  11  is  presented.  Then  the  series  solution  for 
the  plane-wave  scattering  by  a  submerged  fluid-filled  spherical  thin  shell  is  derived.  Next,  the  user’s 
instructions  are  outlined  for  program  SCATSPHERE2.  This  section  includes  discussions  of  the  program’s 
inputs  and  outputs,  as  well  as  its  graphics  capabilities.  Finally,  the  three  examples  presented  are  underwater 
plane-wave  scattering  by  (1)  an  evacuated  steel  spherical  thin  shell;  (2)  a  fluid-filled  steel  spherical  thin  shell; 
and  (3)  a  fluid-filled  aluminum  spherical  thin  shell.  The  first  two  examples  verify  the  updated 
NASHUA/NASTRAN  procedure  while  the  last  example  is  compared  with  (and  shown  to  agree  with)  both 
the  thick  shell  theory  14  and  experimental  results.  15-16 

UNDERWATER  SCATTERING  THEORY 

The  three-dimensional  wave  in  the  fluid  is  modeled  by  the  homogeneous  Helmholtz  equation, 

Ap(x)  +  k2p(\)  =  0,  x  in  R3,  (1) 

where  A  is  the  Laplacian  operator,  p  is  the  pressure,  and  k  is  the  wavenumber.  Here,  the  time-dependent 
term  exp (—iut)  has  been  suppressed.  Using  separation  of  variables  in  spherical  coordinates,  17-18  we 


substitute  p  (r,  9,<J>)=R  (r)6(9)  $(</>),  where  r  is  the  radius,  9  is  the  angle  from  the  z-axis,  and  4>  is  the  angle  in 
the  x-y  plane,  into  Eq.  1  to  obtain  the  spherical  Bessel  equation, 


f(r2  ^)  +  (k2r2-(n2+n))R=0, 
the  associated  Legendre  equation, 

f,  i§-]  +  +")si"2w  - 

and  the  one-dimensional  wave  equation, 

-1  d2<h 

— — - —  =  nr. 

$(d)  dd2 


(2) 


(3) 


(4) 


The  ranges  of  p  and  0  restrict  the  general  solution  to  those  with  integer  eigenvalues,  m  and  n.  The  general 
solution  of  Eq.  2  is  expressed  in  terms  of  j„(kr),  the  spherical  Bessel  function  of  the  first  kind,  and  hn(kr), 
the  spherical  Hankel  function  of  the  first  kind,  both  of  order  n.  Note  that  jn(kr)  is  finite  at  the  origin  and 
unbounded  in  the  infinite  space,  while  hn(kr)  is  unbounded  at  the  origin  and  satisfies  the  Sommerfeld 
radiation  condition  in  the  infinite  space  for  the  assumed  exp (-iut)  time  dependence.  The  general  solution  of 
Eq.  3  is  expressed  in  terms  of  P„(cos9),  the  associated  Legendre  polynomials  of  order  m  and  n.  The  general 
solution  of  Eq.  4  is  the  trivial  one-dimensional  wave  solution.  Note  that  m  is  the  number  of  great  nodal 
circles  intersecting  the  z-axis,  and  n  is  the  number  of  nodal  circles  coaxial  with  the  z-axis.  The  plane-wave 
scattering  by  spherical  obstacles  is  axisymmetric,  and  therefore,  without  loss  of  generality,  di>/d<f>=Q,  m- 0, 
and  ©  solves  the  Legendre  equation. 

Let  the  notation  "J&F"  indicate  a  formula  in  Sound,  Structures,  and  Their  Interaction  by  Junger  and 
Fcit.  11  By  an  addition  theorem,  19(sec,10nl°  1  <l7)  the  free-space  Green’s  function  in  spherical  coordinates  (J&F, 
section  6.10)  is 


(5) 


G(x,x  o)  = 


exp  (i  k  |  x  —  x0 1 ) 
4  7T  |  X  —  X0  I 


=  E  E  +  l)PZ(cos9)PZ(cos00)cos(m(<t>-<t>o))hn(kr)  jn(kr0), 

4~  n=0  m  =0  (»+"»)! 

r  >r0,  x,x0  in  R3, 

where  the  Neumann  factor  (m  is  1  for  m  =  l  and  2  for  m  >1.  By  the  same  addition  theorem,  the  incident 
plane-wave  pressure  (J&F,  section  10.13)  is 


Pair,  6)  =  4> o  E  (2fl  +  l)in  Pn(cos0)]„(kr), 

n  =0 


(6) 


where  4>0  is  the  incident  pressure  magnitude. 


EVACUATED  SPHERICAL  THIN  SHELL 

In  this  section,  the  series  solution  for  the  plane-wave  scattering  by  a  submerged  evacuated  spherical  thin 
shell  is  reviewed.  Assume  that  the  shell  is  thin  enough  and  the  frequency  is  low  enough  that  flexural  stresses 
can  be  ignored  as  compared  to  membrane  stresses.  The  plane-wave  scattering  by  the  submerged  evacuated 
spherical  thin  shell  requires  finding  (1)  the  shell  impedance,  (2)  the  fluid  radiation  impedance  on  the  exterior 
surface  using  the  Green’s  function  for  the  exterior  region,  (3)  the  ’’rigid-body"  scattered  field,  and  (4)  the 
radiated  scattered  field  (due  to  shell  vibration). 

Using  the  equations  for  nontorsional  axisymmetric  motions  (J&F,  sections  7.102  and  7.103),  the  shell 
impedance  for  the  nth  axisymmetric  mode  (J&F,  section  7.121)  is 


Z*  -- 


jjsCp  h  [n2-(nnl)2][n2-(nn2)2] 

n  fl  [n2-(i+/?2)(i/+x„-i)]  ’ 


(7) 


where  ps  is  the  shell  density,  cp  =  V£y[ps(  1-t/2)] ,  E  is  the  Young’s  modulus,  h  is  the  shell  thickness,  a  is  the 
shell  mean  radius,  f 1  =  uia/cp  is  the  dimensionless  frequency,  0  =  h/(a  Vi2),  i'  is  the  Poisson’s  ratio,  and 
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A„  =n  ( n  +1)-  The  quantities  ilnl  and  firt2  are  the  upper  and  lower  shell  resonance  dimensionless  frequencies, 
respectively.  They  are  the  solutions  of  the  characteristic  equation  (J&F,  section  7.114) 

fl4  -fl+3^+An-/i2(l-l,-An2-^An)]n2  +  (A„-2)(1-^)  (8) 

+  /?2  [A„3  —  4  A„2  +  A  (5  —  t/2)  —  2  (1  —  i/2)]  =  0. 


The  exterior  fluid  impedance  (the  fluid  radiation  impedance  at  the  exterior  surface)  is  found  by  first 
finding  the  exterior  region’s  Green’s  function  which  must  satisfy  the  Sommerfeld  radiation  condition  and  the 
Neumann  boundary  condition  The  former  condition  implies  that  the  exterior  region’s  Green’s  function  is  in 
terms  of  h„(kr).  The  latter  condition  and  the  Wronskian  relation  19(secIlons81411ana11414)  implies  that  the 
Green’s  function  for  the  exterior  region  specialized  '  i  the  spherical  surface  (J&F,  section  6.15)  is 


1  “  «  (n—mV  hn(kr) 

G(r,0,d\a,eo,i>o)  =  — f-r  E  E  <m  (2n  +  l)cos(m(4-4o))PZ{cosO)P?(cos00)  -~TT>  (9) 

4ffka-„=n„,=n  («+«)!  K  (ka) 


n  =0  m  =0 


(n  +m)l 


r>a. 


The  Green’s  identity  for  the  exterior  region  and  the  problem’s  axisymmetry  imply  that  the  fluid  pressure 
(J&F,  section  6.16  and  6.19)  is 


p(r,e,<t> )  =  —  J '  G{r,6,<t>\a,0o,<l>o)w'{6o,<t>o)sin{eo)dS 
ioj  Js 

00  ,  hn{kr) 

=  ipfCf  E  P”(cos0)TG7rT’  T>0’ 

n=0  \kQ) 


(10) 


where  is  the  fluid  density,  w'  is  the  shell  surface  normal  velocity,  and  Cf  is  the  fluid  sound  speed.  The 
exterior  fluid  impedance  of  axisymmetric  mode  n  (J&F,  section  6.29)  is 


^•n 


p„(a,9,<t>)  _  .  hn(ka ) 
wn'(a,6,<t>)  1  PfCf  hn'(ka) ' 


(11) 


The  "rigid-body"  scattered  pressure  field  (J&F,  section  10.16)  is  found  by  using  the  Neumann  boundary 
condition  and  the  incident  plane-wave  pressure.  Using  the  asymptotic  expansion  for  the  Bessel  function,  the 


far-field  "rigid-body"  scattered  pressure  field  U&F,  section  10.17)  is 


ie‘kr<b o  00  jn'(ka'' 

pAr,0,i>)  =  - — ZT~  £  (2n+l)Pn(cosO)  r»a. 


kr 


n  =0 


hn\ka) 


(12) 


From  (J&F,  section  11.25)  to  (J&F.  sectio  i  11.28),  the  shell  velocity  response  to  the  incident  and 
"rigid-body"  scattered  pressure  is 


„  „  .X  ™PnPn(cOS0) 

w  (r,0,o)  =  -  £  (  , 


(13) 


<1>0  in+1  (2/1  -rl)Pn(cosO) 
k2a2  „?0  (Zn  +  z„)hn\ka) 


where  p„=(Zn  +  z„)  wn'  is  used.  From  Eqs.  10  and  13  and  the  asymptotic  expansion  for  the  Bessel  function, 
the  far-field  radiated  pressure  due  to  the  shell  motion  (J&F,  section  11.29)  is 


,  „  i  Pfcfe  k  <I>o  ^  {2n+l)Pn(cos6) 

pr{r,0,6)  =  —  — -  £  — — —— — -  -'T  >  r»a. 

kr  „=0(Zn+zn)[kahn(ka)Y 


(14) 


The  total  far-field  scattered  pressure  (J&F,  section  11.30),  which  is  the  sum  of  the  "rigid-body"  scattered 
pressure  and  the  radiated  pressure,  is 


ie‘kr<t> o  00  (2n  +l)Pn(cos6) 
Pse{r,0,<S>)  =  - — — —  £ 


kr 


n  »0 


hn'(ka) 


j'n(ka)  - 


Pfcf 


(Zn+zn)(ka)2hn'(ka) 


,  r»a. 


(15) 


For  more  in-depth  theoretical  presentation,  see  Sound,  Struiiures,  and  Their  Interaction,  by  Junger  and 
Feit.  11 
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FLUID-FILLED  SPHERICAL  THIN  SHELL 

In  this  section,  the  analytical  plane-wave  scattering  solution  for  the  submerged  fluid-filled  spherical  thin 
shell  is  derived  by  extending  the  theory  from  the  previous  section.  The  extended  solution  requires  finding  (1) 
the  interior  fluid  impedance  using  the  Green’s  function  for  the  interior  region,  and  (2)  the  radiated  pressure 
field  due  to  shell  vibration. 

The  Green’s  function  for  the  interior  region  is  expressed  in  terms  of  j„{kr)  because  it  must  be  finite 
inside  the  sphere.  Using  the  Neumann  boundary  condition  to  find  the  unknown  coefficient,  the  Green’s 
function  for  the  interior  region  is 


G{r,0,6  |fl,0o.bo)  = 


1  oo  OO  (n—m  V 

— r  E  E  -T7  (2»+l)cos(m  (b- 

4  "  &  a*  n=0m=0  («+"»)! 


■<l>o))Pn(COSO)P”(cOS00) 


-jn(kr) 

jn'(ka) 


,  (16) 


r  <a. 


'1'he  Green’s  identity  for  the  interior  region  and  the  problem’s  axisymmetry  imply 
p{r,0,<t>)  =  —  J  G(r,e,<p\a,Oo,^o)w'(Oo,<j>o)sin(0o)dS 

l  U) 


(17) 


OO 


=  i  PgCg  E  wn'  Pn(cos0) 

n  =0 


-jn(kr) 
jn'(ka)  ’ 


r  <a, 


where  pg  is  the  interior  fluid  density,  and  cg  is  the  interior  fluid  sound  speed.  Thus,  the  interior  fluid 
impedance  of  mode  n  is 


C/I 


Pn 


=  t 


Pg  Cg 


-]n(ka) 

in(ka)  ' 


(18) 


The  only  difference  between  computing  the  evacuated  and  fluid-filled  cases  is  that  the  interior  fluid 
impedance  fn  needs  to  be  added  to  the  shell  modal  impedance  Z„.  Thus,  the  far-field  radiated  pressure  is 
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pr(r,6,<f>)  =  - 


ipfCfe‘kr4> o  00  (2n  +  l)Pn(cos6) 


kr 


^  (Zn+z„+$n)[kahn'(ka)]2  ’ 


(19) 


n  =0 


The  total  far-field  scattered  pressure  (which  is  the  sum  of  the  "rigid-body”  scattered  pressure  and  the 
radiated  pressure)  is 


P*{r,0,<j>)  = 


ie‘kr<t> o  00  (2n+l)Pn{cos6) 
kr  n5o  hn'(ka) 


I  i’nika) 


_ Pf£f _ 

(Zn+zn  +  $n)(ka)2  hn\ka ) 


r  »  a. 


(20) 


PROGRAM  IMPLEMENTATION 

Program  SCATSPHERE2  computes  the  spherical  Bessel  function  values  after  using  IMSL  library  calls 
for  the  cylindrical  Bessel  function  values.  19(sectiomo.i.i)  -p^e  Legencire  polynomials  are  computed  using  a 
recursion  relation.  20  The  infinite  series  is  approximated  by  a  finite  number  of  terms.  More  terms  are  added 
until  either  the  relative  error  criterion  is  met  or  the  maximum  number  of  terms  allowed  is  used. 

USER’S  INSTRUCTIONS 

SCATSPHERE2  input  data  consist  of  11  record  types  for  25  parameters  which  are  illustrated  in  Table 
1.  All  variables  are  read  in  as  unformatted  (i.e.,  separated  by  commas).  Records  7a  and  7b  are  read  in  only 
if  the  field  pressure  calculations  are  done,  i.e.,  NFLAG  =  2  or  3.  Record  7b  consists  of  an  array  of  radii 
whose  dimension  depends  on  NRADII.  After  NINT  is  recorded, the  next  record,  Record  11,  is  read  NINT 
times.  The  dimensionless  frequency  in  Table  1  is  defined  as  fi  =  k  a  Cf/cp.  The  complex  Young’s  modulus  is 
E(l  +  i  rj). 

The  series  is  computed  for  as  many  terms  as  needed  until  the  relative  error  criterion  is  satisfied  or  the 
maximum  number  of  terms  is  used.  Recommended  values  are  NTERMS  =  40  and  ERRCRT  =  1.0  x  10-6.  If 
erratic  behavior  occurs,  then  increase  NTERMS  and  decrease  ERRCRT.  Overall,  the  computation  is  quite 
inexpensive  even  for  small  ERRCRT.  The  maximum  value  for  NTERMS  is  201,  and  the  maximum  number  of 
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Table  1.  SCATSPHERE2  input  data. 


Description 

1 

Title  Card  (up  to  75  characters) 

o 

$0 

Incident  Pressure  Magnitude 

3a 

a 

Shell  Radius 

3b 

h 

Shell  Thickness 

3c 

Ps 

Shell  Density 

3d 

V 

Shell  Poisson’s  Ratio 

3e 

E 

Shell  Young’s  Modulus 

3f 

V 

Shell  Damping  (or  Loss)  Factor 

4a 

Pf 

Exterior  Fluid  Density 

4b 

cf 

Exterior  Fluid  Sound  Speed 

5a 

Ps 

Interior  Fluid  Density  (=  0  for  evacuated  shell) 

5b 

r 

Lg 

Interior  Fluid  Sound  Speed  (=  0  for  evacuated  shell) 

6a 

IPRINT 

Print  Control  Flag  (0-5) 

0  =  Minimal  Printing 

1  =  Normal  Printing 

2  =  Above  Normal  Printing 

3  =  Low  Debug  Printing 

4  =  Debug  Printing 

5  -  Excess  Printing 

6b 

NFLAG 

Calculation  Type  Flag  (1-5) 

1  =  Surface  Pressure  Calculation  Only 

2  =  Surface  and  Field  Pressure  Calculation 

3  =  Field  Pressure  Calculation  Only 

4  =  Surface  Pressure  and  Radial  Velocity  Calculation 

5  =  Surface  Radial  Velocity  Calculation  Only 
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Table  1  (Continued) 


Record  No. 

Parameter 

Description 

7a 

NRADII 

Number  of  Observation  Radii:  If  NFLAG  =  2  or  3 

7b 

r 

Radii  Array:  1  to  NRADII:  If  NFLAG  =  2  or  3 

8a 

NTERMS 

Maximum  Number  of  Series  Terms 

8b 

ERRCRT 

Relative  Error  Criterion 

9a 

h 

Starting  Angle  from  z-Axis 

9b 

A  6 

Delta  Angle  from  z-Axis  (>  0) 

9c 

02 

Ending  Angle  from  z-Axis  (>  0X) 

10 

NINT 

Number  of  Frequency  Runs 

11a 

fii 

Starting  Dimensionless  Frequency:  1  to  NINT 

lib 

An 

Delta  Dimensionless  Frequency:  1  to  NINT 

11c 

n2 

Ending  Dimensionless  Frequency:  1  to  NINT 

samples,  which  is  the  number  of  ka' s  times  the  number  of  angles,  is  2001. 

Program  SCATSPHERE2  produces  two  output  files,  tapes  11  and  13.  Tape  11  contains  the  input  data 
and  a  table  for  the  nondimensional  pressure  \pr  /p0a  | .  Tape  13  contains  an  input  X-Y  PLOT  file  which  is 
processed  by  the  LOGOS  graphics  program.  21  Tape  12  is  not  currently  used,  but  the  user  could  restore  this 
tape  for  input  to  F.M.  Henderson’s  program  PLOTTER,  which  creates  polar  plots  and  uses  the  software 
DISSPLA.  22  For  more  information  on  using  program  PLOTTER,  see  F.M.  Henderson’s  report  on 
RADSPHERE.  12 

The  setup  for  a  typical  CRAY  COS  job  is  as  follows: 

JOB,JN=SCAT,T=60,US=BUDGET.  RUN  SCATSPHERE 
ACCOUNT,  AC=xxxxxxxxxx,US=yyyy,UPW=zzzz. 

CFT77,ON=AXM. 

ACCESS, DN=IMSL,OWN=PUBLIC. 

SEGLDR,CMD=’LIB=IMSL’,GO. 
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DISPOSE, DN=TAPE11,... 
DISPOSE, DN=TAPE12, . . . 
DISPOSE, DN=TAPE13, . .. 

/EOF 

PROGRAM  SCATSPHERE2 


END 

/EOF 

Forward  Scattering  for  a  Fluid-Filled  Steel  Spherical  Shell 

20. 

5. . 0.15.7669.. 0.3. 2.07E11, 0.01 

1000.,  1524. 

1000..  1524. 

1,3 

1,100. 

80,1.E-S 

0.,30.,0. 

1 

0.005,0.005,5.0 

Note  that  the  code  is  in  FORTRAN-77  and  IMSL  library  is  used  for  the  Bessel  function  calculations. 


NUMERICAL  EXAMPLES 

In  this  section,  we  present  three  examples:  underwater  plane-wave  scattering  by  (1)  an  evacuated  steel 
spherical  thin  shell,  (2)  a  fluid-filled  steel  spherical  thin  shell,  and  (3)  a  fluid-filled  aluminum  spherical  thin 
shell.  The  first  two  examples  verify  the  updated  NASHUA  procedure,  while  the  last  example  is  compared  to 
and  agrees  with  the  numerical  implementation  of  the  thick  shell  theory  14  and  experimental  results.  15-16  All 
the  examples  are  for  a  thin  shell;  that  is,  the  shell  thickness  is  no  more  than  5%  of  the  shell  radius. 


EVACUATED  STEEL  SPHERICAL  THIN  SHELL 

This  example  was  used  to  verify  the  updated  NASHUA  procedure  and  to  check  program 
SCATSPHERE2.  Figure  1  illustrates  the  time-harmonic  plane-wave  scattering  problem  for  a  submerged 
evacuated  steel  spherical  thin  shell.  Table  2  presents  the  shell  and  fluid  properties  used.  Note  that  the 
interior  fluid  properties  pg  and  cg  are  set  to  zero. 
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Fig.  1.  Plane-wave  scattering  from  a  submerged  elastic  spherical  shell. 


Table  2.  Spherical  thin  steel  shell  and  fluid  properties. 


Description 

Parameter 

Measurement 

Shell  Radius 

a 

5.0  m 

Shell  Thickness 

h 

0.15  m 

Shell  Density 

Ps 

7669.0  kg/m 3 

Shell  Poisson  Ratio 

V 

0.3 

Shell  Young  Modulus 

E 

2.07x10"  M/m2 

Shell  Loss  Factor 

V 

0.0 

Exterior  Fluid  Density 

Pf 

1000  kg/m3 

Exterior  Fluid  Sound  Speed 

cf 

1524  m  /s 

Interior  Fluid  Density 

Ps 

0  kg/m 3 

Interior  Fluid  Sound  Speed 

cs 

0  m  /s 

For  the  NASHUA  finite  element  model,  40  axisymmetric  conical  shell  elements  spanning  the  180 
degrees  between  the  two  spherical  poles  were  used,  making  a  total  of  205  structural  and  41  fluid  degrees  of 
freedom  (DOF)  for  the  exterior  fluid.  The  nondimensional  frequency  increment  was  about  ka  =0.02  in  the 
range  0<ka  <5. 
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For  SCATSPHERE2:  the  maximum  number  of  series  terms  NTERMS  =  40,  the  relative  error  criterion 
ERRCRT  =  1.0x10^,  NFL  AG  =  3  for  the  field  pressure  calculation  only,  N1NT  =  1  for  one  frequency 
sweep,  and  NRADII  =  1  for  one  radial  observation  per  frequency  point.  The  nondimensional  frequency 
range  0 <ka  <5  was  swept  using  a  frequency  increment  of  ka  =0.005. 

Figures  2  and  3  show  the  forward  and  back-scattered  pressure  fields  vs.  the  dimensionless  frequency. 
These  figures  show  that  NASHUA  and  SCATSPHERE2  agree  well. 

FLUID-FILLED  STEEL  SPHERICAL  THIN  SHELL 

This  example  was  used  to  verify  the  updated  NASHUA  procedure.  The  properties  are  the  same  as  in 
the  previous  example  except  that  the  spherical  thin  shell  is  filled  with  the  same  interior  fluid  as  the  exterior 
fluid. 

For  the  NASHUA  finite  element  model,  40  axisymmetric  conical  shell  elements  spanning  the  180 
degrees  between  the  two  spherical  poles  were  used,  making  a  total  of  205  structure  DOF  and  41  fluid  DOF 
for  each  fluid  region.  The  nondimensional  frequency  increment  was  about  Am  =0.05  in  the  nondimensional 
frequency  range  0 <ka  <5. 

For  SCATSPHERE2,  the  maximum  number  of  series  terms  NTERMS  =  80,  the  relative  error  criterion 
ERRCRT  =  1.0x10”^,  NFLAG  =  3  for  the  field  pressure  calculation  only,  NINT  =  1  for  one  frequency 
sweep,  and  NRADII  =  1  for  one  radial  observation  per  point.  Note  that  NTERMS  and  ERRCRT  are 
different  from  those  used  in  the  evacuated  case.  These  changes  were  made  to  avoid  convergence  problems 
due  to  the  Bessel  function  computations.  The  nondimensional  frequency  increment  was  ka  =0.005  in  the 
nondimensional  frequency  range  0 <ka  <5. 

Figures  4  and  5  plot  the  frequency  response  of  the  forward  and  back-scattered  nondimensional  pressure 
\pr  /p0a  | .  Both  figures  show  very  good  agreement  between  the  two  approaches. 
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Fig. 


K  A 

2.  Forward-scattered  field  by  a  submerged  evacuated  steel  spherical  shell. 


K  A 


Fig.  3.  Back-scattered  field  by  a  submerged  evacuated  steel  spherical  shell. 


K  A 

Fig.  4.  Forward-scattered  field  by  a  submerged  fluid-filled  steel  spherical  shell. 


K  A 

Fig.  5.  Back-scattered  field  by  a  submerged  fluid-filled  steel  spherical  shell. 
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Table  3.  Spherical  thin  aluminum  shell  and  fluid  properties. 


Description 

Parameter 

Measurement 

Shell  Radius 

a 

4.875  m 

Shell  Thickness 

h 

0.25  m 

Shell  Density 

Ps 

2700.0  kg/m 3 

Shell  Poisson  Ratio 

V 

0.355 

Shell  Young  Modulus 

E 

0.675  xl0u  N/m 2 

Shell  Loss  Factor 

V 

0.00 

Exterior  Fluid  Density 

Pf 

1000  kg/m 3 

Exterior  Fluid  Sound  Speed 

cf 

1410  m/s 

Interior  Fluid  Density 

Pg 

1000  kg  /m  3 

Interior  Fluid  Sound  Speed 

cg 

1410  m/s 

K  A 

Fig.  6.  Back-scattered  field  by  a  submerged  fluid-filled  aluminum  spherical  shell. 
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FLUID-FILLED  ALUMINUM  SPHERICAL  THIN  SHELL 

This  example  was  chosen  from  earlier  papers  14-15  as  another  example  to  verify  program 
SCATSPHERE2.  Table  3  presents  the  shell  and  fluid  properties  used.  Figure  6  gives  the  back-scattered 
nondimensional  absolute  pressure  for  nondimensional  frequency  0 <ka  <20.  Close  comparison  of  Fig.  6  with 
Fig.  4b  in  Hickling’s  paper  14  and  Fig.  lb  in  Diercks’s  and  Hickling’s  paper,  15  shows  good  agreement 
between  SCATSPHERE2  computations  and  results  obtained  by  Diercks  and  Hickling. 

SUMMARY 

Program  SCATSPHERE2  solves  the  time-harmonic  plane-wave  scattering  by  a  submerged  elastic 
spherical  thin  shell  that  is  evacuated  or  fluid-filled.  This  code  was  verified  by  the  updated  NASHUA 
procedure,  24  the  numerical  implementation  of  the  thick  shell  theory,  14  and  experimental  results.  15  This 
program  allows  the  user  to  examine  the  surface  velocity  at  any  angle,  and  the  scattered  pressure  field  at  any 
radius  and  angle.  Program  SCATSPHERE2  also  allows  the  user  to  create  polar  plots,  tables,  and  X-Y  plots. 
In  the  future,  program  SCATSPHERE2  could  be  modified  for  other  axisymmetric  problems  such  as  including 
additional  interior  mediums  (solids  and/or  fluids)  or  spherical-wave  scattering. 
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1 

1941 

J.  Su 

1 

1282  G.C.  Everstine 

1 

1941 

R.  Vasudevan 
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CENTER  DISTRIBUTION  (Continued) 


Copies 

Code 

Name 

1 

1941 

W.H.  Vogel 

1 

1942 

1 

1942 

Y.-F.Hwang 

1 

1943 

1 

1944 

1 

1945 

1 

1945 

K.  Jones 

1 

27 

1 

2704 

1 

2704.1 

1 

274 

1 

2741 

1 

2740 

D.E.  Goldsmith 

1 

2742 

1 

2742 

H.C.  Neilson 

1 

2743 

1 

2744 

1 

2749 

1 

28 

1 

342.1 

(TIC) 

1 

342.2 

(TIC) 

10 

3432 

Reports  Control 

1 

35 

S.E.  Willner 

1 

93 
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